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The normal spinel oxide MgCr2C>4 is cubic at room temperature while the normal spinel CuC^CU 
is tetragonal as a consequence of the Jahn-Teller activity of Cu 2+ on the tetrahedral sites. Despite 
different end-member structures, a complete solid solution of Mgi-^CuzC^CU can be prepared with 
compounds of composition x = 0.43 displaying a first-order phase transition at room temperature. 
Reverse Monte Carlo analysis of total neutron scattering on data acquired between 300 K and 15 K 
on samples with x = 0.10, 0.20, and 0.43 provides unbiased local and average structure descriptions 
of the samples, including an understanding of the transition from local Jahn-Teller distortions 
in the cubic phase to cooperative distortions that result in a tetragonal structure. Distributions 
of continuous symmetry measures help to understand and distinguish distorted and undistorted 
coordination around the tetrahedral site in the solid solutions. Magnetic exchange bias is observed 
in field-cooled hysteresis loops of samples with dilute Cu 2+ concentration and in samples with 
tetragonal-cubic phase coexistence around 300 K. 



PACS numbers: 71.70.Ej, 61.05.fm, 75.50.Tt 

INTRODUCTION 

The propensity of octahedral Cu 2+ ions in oxide struc- 
tures to display Jahn-Teller (JT) distortions is intimately 
linked to magnetism and superconductivity in systems 
derived from La2Cu04. [HE] While less common, tetra- 
hedral Cu 2+ on the A site of oxide spinels can also display 
JT activity. This distortion lowers symmetry by com- 
pressing the tetrahedron and thereby breaking the degen- 
eracy of the partially-occupied ti energy levels. [2l|4] The 
crystal field splittings of the ideal and distorted tetrahe- 
dra are shown schematically in Fig. [I] 

Jahn-Teller distortions themselves are an intriguing 
theme in functional materials because they enable inter- 
play between electronic and structural degrees of free- 
dom. They have been most widely studied in the man- 
ganites, often derivatives of perovskite LaMnC-3. [SHI] In 
these compounds, Mn 3+ has four 3d electrons with a 
singly-occupied pair of e g states in an octahedral crystal 
field. It is well established that elongation of the octahe- 
dron breaks the degeneracy and lowers the energy of the 
system. [5] [9] The percolative nature of orbital ordering 
arising from cooperative JT distortion is believed to play 
a central role in colossal magnetoresistive behavior. [10] 

In spinels, collective JT distortions on the A or B sites 
result in a reduction in symmetry from cubic Fd3m to 
tetragonal IAi/amd upon orbital ordering at the Jahn- 
Teller transition temperature Tjt- When only a fraction 
of occupied sites are JT-active, cation clustering can lead 
to endotaxial coexistence of tetragonal (distorted) and 
cubic phases, [TTHE] with strain-driven checkerboard pat- 
terns first studied in phase-separated CoPt alloys. [T4T 
[T6] These Mn 3+ -driven JT distortions are a product of 
unpaired 3d electrons, so self-assembled nanostructured 
magnetic films are under development. [17] 
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FIG. 1: (Color online) Crystal-field splitting of Cu 2+ with a 
3d 9 electron configuration in an ideal tetrahedron (left) results 
in degeneracy in the t^ orbitals. Jahn-Teller compression of 
the tetrahedron lifts the degeneracy (right), and results in a 
singly-occupied xy orbital. 



Few studies have examined the precise JT tendency of 
Cu 2+ on the spinel A site. The effect of Cu 2+ occupancy 
on Tjt and the electronic or magnetic properties remains 
sparesly investigated. [151120] We expect there should 
be key differences between JT activity on the A and B 
sites of spinel. While BOq octahedra are edge-sharing 
and form a pyrochlore sublattice, the AO4 tetrahedra are 
isolated from each other in a diamond sublattice. The 
increased distance between A cations should hinder their 
cooperative behavior. 

In this work, the structural effects of Cu 2+ concen- 
tration x in the spinel solid solution Mgi-zCuzC^C^ 
are investigated by real- and reciprocal space structural 
probes utilizing total neutron scattering. Discrepancies 
between the average and local structural are of partic- 
ular interest. Specifically, we probe whether the coor- 
dinations of JT-active Cu 2+ and JT-inactive Mg 2+ dif- 
fer. Traditional Bragg diffraction analysis fails to resolve 
these differences because atoms on the same crystallo- 
graphic site (here Cu and Mg on the spinel A site) are 
required to have identical surroundings. Electron para- 
magnetic resonance [2TJ and x-ray absorption [55] offer 
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a confirmation that dilute Jahn- Teller cations create lo- 
cal distortions, but these spectroscopic techniques do not 
yield any detailed structural information. We utilize the 
pair distribution function (PDF) because it provides a 
real-space description of the structure with distinct Cu- 
O and Mg-0 distances. 

Our previous PDF study of the spinel CuMn 2 04 
encountered many of the complications that make 
Mgi-zCu^C^CU a difficult crystal structure to de- 
scribe. [23 In that study, Cu and Mn are present with 
mixed valence on both A and B positions, and the sur- 
rounding oxygen polyhcdra arc cation-dependent. We 
found that CUO4 tetrahedra are significantly more dis- 
torted (as judged by bond angles) than MnC>4 tetra- 
hedra, and Cu avoids the tendency for JT distortion 
by disproportionating to Cu + /Cu 3+ . The situation for 
Mgi-zCu^C^C^ should be less complex: no valence 
mixing is present, and Cu/Mg substitution is confined 
to the A site. The effect of central cation on MO4 dis- 
tortion is more isolated. 

We employ large box modeling via reverse Monte Carlo 
simulations as method of retrieving possible signatures of 
cation-dependent coordinations from the PDF. [2H [S3] 
Supercells with thousands of atoms can be routinely 
simulated with modest computational requirements, and 
the large sample size provides element-specific informa- 
tion due to the presence of many discrete atoms of each 
type. A crucial aspect of RMC simulations is deter- 
mining straightforward metrics that describe how local 
crystalline structure differs from the average. The ten- 
dency for distortion of individual polyhedra has been 
characterized by analyzing bond lengths, [53] bond an- 
gles, [H W\ and geometric analysis. gl US] 

In this study, tetrahedral JT distortion is gauged us- 
ing continuous symmetry measures (CSM). [30] [3T] The 
particular strength of the CSM method is its ability to 
compare the symmetry of imperfect polyhedra, regardless 
of their size or orientation in space. |321 133j Extraction 
of CSM information from RMC simulation was recently 
performed as a test of structural rigidity in Bi 2 Ti 2 0y [34] 
and is employed here to compare the symmetry of CUO4, 
Mg04, and AA4 tetrahedra, thereby describing the pref- 
erence for JT activity as a function of x and temperature. 

METHODS 

Powders of Mgi-zCuzC^C^ compounds were 
prepared by dissolving stoichiometric amounts 
of Cu(N0 3 ) 2 -2.5H 2 0, Mg(N0 3 ) 2 -6H 2 0, and 
Cr(N03)3-9H 2 in water, followed by boiling to 
evaporate the solvent until a brown mass was formed, 
which was then ground and calcined in air at between 
700°C or 1000°C for 10 hours, then cooled at 10°C/min. 
Laboratory X-ray diffraction patterns were acquired us- 
ing Cu-K a radiation on a Philips X'Pert diffractometer 



at room temperature and a Bruker D8 diffractometer 
with an Anton Parr high-temperature stage. Magnetic 
properties were measured using a Quantum Design 
MPMS 5XL SQUID magnetometer. Time-of-flight 
(TOF) neutron scattering was performed on the NPDF 
instrument at Los Alamos National Laboratory. Rietveld 
refinements were performed using the XND code [35] 
for X-ray data and the GSAS-EXPGUI suite[36] for 
the TOF data. The PDF was extracted using the 
PDFGetN program [57] with Q max = 35 A -1 and least- 
squares refinement of the PDF was performed using the 
PDFgui frontend for PDFFit2. 38j Crystal structures 
were visualized with AtomEye[39] and VESTA. [40] 

Reverse Monte Carlo simulations were run using RM- 
CProfile version 6 [H] on 7 x 7 x 7 cubic or 10 x 10 x 7 
tetragonal spinel supercells with 19208 or 19600 atoms, 
respectively. A hard-sphere repulsion was applied to pre- 
vent M-0 bond distances shorter than the first peak of 
the PDF, but no clustering was observed at the cutoff 
distances. No preference for cation clustering was found, 
so configurations with randomized Cu and Mg occupancy 
were used. Simulations were performed as serial jobs on 
the HP Opteron QSR cluster at the California NanoSys- 
tems Institute. 

Bond valence sums (BVS) were extracted from atoms 
in the superccll in the same manner described in our pre- 
vious work on CuMn9Q^.[23] using the Ro values of Brese 
and Q'Keeffe.[4"2"] CSM for A0 4 tetrahedra were calcu- 
lated using a distance measure program provided by M. 
Pinsky and D. Avnir. 

RESULTS AND DISCUSSION 
Average structure via reciprocal-space analysis 

The compounds MgCr 2 04 and CuCr 2 04 both belong 
to the AB2O4 spinel family of structures with the A 
cations, Mg 2+ and Cu 2+ , tetrahedrally coordinated by 
oxygen, while Cr 3+ lies on the octahedral B site. The 
[Ar]3<i 3 electron configuration of Cr 3+ is very stable be- 
cause each of the t 2g energy levels is singly occupied, so 
there is no tendency of site mixing or mixed valence. [43] 
Varying x in Mgi_ 2; Cu 2 ;Cr 2 04 therefore does not disturb 
the B sublattice composition, but changes in the inter- 
penetrating A sublattice may lead to chemical pressure 
which will influence its size and shape. 

On the tetrahedral A site, Mg 2+ and Cu 2+ have ef- 
fectively identical ionic radii. Both are 0.57 A as given 
by Shannon, [44] but their electron configurations are dis- 
tinctly different. Mg 2+ has the [Ne] configuration and no 
d electrons. Cu 2+ has [Ar]3<i 9 and only two of the three 
ti energy levels are fully occupied in tetrahedral coor- 
dination (Fig.[l]). This de generacy causes a JT distor- 
tion, manifested by a flattening of the tetrahedron. Bond 
lengths are preserved, but bond angles are no longer 
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FIG. 2: (Color online) To scale, the pseudotetragonal cells of 
(a) cubic MgCr2C>4 and (b) Jahn-Teller distorted tetragonal 
CuCr2 04 are viewed along the a axis of the IAi/amd cell. 
Contraction in the c direction is evident due to JT distortion 
in CuCr204. Cr06 octahedra are light blue, while Mg04 and 
CUO4 tetrahedra are dark blue and red, respectively. 
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FIG. 3: (Color online) Lattice parameters of the 
Mgi_ a: Cu :E Cr2 04 solid solution obtained by Rietveld refine- 
ment of laboratory X-ray powder diffraction at room temper- 
ature. Coexistence of the cubic and tetragonal phases occurs 
for 0.43 < x < 0.47. V 1 ^ 3 is shown for tetragonal phases. 
Error bars are smaller than the symbols for all points. 



equivalent at 109.5°. [45H47] The contrasting behavior of 
Mg 2+ and Cu 2+ in tetrahedral coordination is evident 
when the MgCr 2 04 and CuCr 2 04 structures are com- 
pared in Fig. [2] MgCr 2 04 forms in the cubic space group 
Fdim with ideal Mg04 tetrahedra. [15] CuCr 2 04 under- 
goes flattening in the c direction and forms in the tetra- 
hedral space group IAx/amd, which is the same space 
group as JT distorted Mn 3 04 (and other AMJI2O4) or 
NiCr9Q4.|45l FI9"] The unit cells are shown to scale in 
Fig. [2] to highlight their difference in dimensions. 

Alloying x from to 1 in Mgi_ a; Cu a; Cr 2 04 produces a 
transition from a cubic to tetragonal spinel. The lattice 
parameters obtained from Rietveld refinements to room- 
temperature XRD patterns shown in Fig. [3] reveal that 



800 



600 



400 



200 







A' 



Cubic 



h A'' 

0/ 1/ 



/ / 
/ / 
/ / 
/ / 
/ / 



□ □ 

/ 1/ 

/ & 

n ' ' 

/ / 
/ / 

1 1 



Tetragonal 



600 



400 



200 " 



-200 



MgCr 2 4 0.2 0.4 0.6 0.8 CuCr 2 4 
xin CuM gi . x Cr 2 4 

FIG. 4: (Color online) Phase diagram of the Mgi_ a; Cu :E Cr2 04 
system as determined by Rietveld refinement. Points denote 
neutron refinements to cubic (□) and tetragonal (x) phases. 
Bars represent coexistence regions from high-temperature X- 
ray diffraction. Letters (a,b,c) correspond to the Rietveld 
refinements shown in Fig. [5] 



for x < 0.43 a cubic spinel is formed with a gradually de- 
creasing lattice parameter a. When x = 0.43 the tetrag- 
onal phase appears with c/a = 0.975 and a small region 
of coexistence persists for 0.43 < x < 0.47, above which 
the cubic phase disappears. The compound becomes in- 
creasingly tetragonal as the end member CuCr 2 04 is ap- 
proached, with c/a — 0.912 when x = 1. The pseudocu- 
bic cell volume contracts from 579.4 A 3 for MgCr 2 04 to 
566.1 A 3 for CuCr 2 4 , which is a 2.3% decrease. 

We extend our Rietveld analysis using high- 
temperature XRD and low-temperature TOF neutron 
scattering to produce an approximate phase diagram of 
the pseudobinary system MgCr 2 04-CuCr 2 04, shown in 
Fig. [4] The Tjt between cubic and tetragonal spinels has 
a nearly linear relationship on x. There is some phase co- 
existence determined from HTXRD denoted by the bars 
on the graph. We use the end-member CuCr 2 04 transi- 
tion of T JT = 590°C from the literature. [Ml [50] The JT 
transition temperature steadily decreases with decreasing 
x so that the tetragonal phase occurs for x = 0.20 but 
not for x = 0.10 at T — 15 K (the lowest temperature 
measured). Rietveld refinements to neutron scattering 
data are shown in Fig. [5] for three different points on the 
phase diagram, representing (a) cubic, (b) mixed, and (c) 
tetragonal phases at x = 0.20 and T = 300 K, x = 0.43 
and T = 300 K, and x = 0.20 and T = 15 K, respectively. 
In all cases, the fits are excellent. There is, however, some 
unfit intensity between split tetragonal peaks in (b) and 
(c), indicative of a more complex crystal structure than 
the two-phase Rietveld model would suggest. 

Not shown on our phase diagram is the tetragonal dis- 
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FIG. 5: (Color online) Time-of-flight neutron scattering Ri- 
etveld refinements of Mgi-zCu^C^CU with (a) x — 0.20 at 
300 K [cubic], (b) x = 0.43 at 300 K [cubic and tetragonal], 
and (c) x = 0.20 at 15 K [tetragonal]. 



tortion around 12 K in MgCr^O^. [5TH53"] A transition 
from FdZm to IAi/amd occurs to remove the geomet- 
ric frustration of Cr 3+ spins arranged in a pyrochlorc 
sublattice. Addition of magnetic Cu 2+ cations in the A 
sites should relieve this frustration due to strong A-B in- 
teractions, [T51 [M] so we do not expect the spin-driven 
distortion to play a role for x > 0. The sample with x — 
0.10 is cubic at 15 K. 

The presence of a tetragonal phase for low Cu 2+ con- 
tent is surprising. At x = 0.20, for instance, only one in 
five A sites has a JT active cation. The A cations are ar- 
ranged in a diamond sublattice and each has four nearest 
A neighbors. The tetrahedra do not share edges or cor- 
ners with each other, with the shortest exchange pathway 
being A— O-Cr— O— A Given a random cation distribu- 
tion, the probability of one Cu 2+ having all JT-inactive 
Mg 2+ neighbors is (f) 4 = 41.0%. The probability of 
having only one Cu 2+ neighbor is 4(|)(|) 3 = 41.0%, 
and the probability of having two Cu 2+ neighbors falls 
to 6(i) 2 (|) 2 = 15.4%. Thus 82.0% of Cu 2+ cations have 
zero or only one JT-active nearest neighbor, but they 
still produce orbital ordering with long-range periodic- 
ity. The trend of Tjt versus x is roughly linear in Fig. [4] 
with no apparent jump at the percolation threshold of 
the diamond-type A sublattice. [551 [55] 

The critical concentration of Cu 2+ needed to drive a 
cooperative JT distortion in Mgi-zCiizCh^C^ at 300 K 
is x = 0.43. This fraction increases with Asite cation 
radius in ACT2O4 spinels: Zni^CuzC^O^ (r\z„ = 0.60 
A) is reported to have x = 0.47 [57] and x = 0.58, [T5] 
while Cdi_zCuzCr 2 4 ( r cd = 0.78 A) has x = 0.64. [2D] 
This could be due to increased distance between Asite 



cations, or a loosening of the structure (thus weakening 
of strain field produced by a JT distortion). 

The critical concentration in Zni-.rNi^C^O^, [58] 
where Ni 2+ drives JT disortion, is around x ~ 1 at room 
temperature. For Cr 3+ on the spinel B site, less Cu 2+ 
is needed to drive a cooperative distortion than Ni 2+ . 
The 3d 8 configuration of Ni 2+ has only one unpaired ti 
electron, rather than the two of Cu 2+ . The result is a 
smaller energy gain after breaking degeneracy and elon- 
gation (rather than contraction) of the c axis. [J] 

Comparison with JT tendency of Mn 3+ on the spinel 
B site is less direct. For example, the solid solution 
ZnfFei-^Mn^O^ has a critical concentration of about x 
= 0.3, [5j] while Mn[Cri_ :c Mn K ] 2 04 has x = OA|60J This 
would seem to indicate a stronger JT tendency, in part 
due to closer B-B distances and edge sharing between 
octahedra. However, for Zn^Ge^^ICoi-zMn^O^, 
Wickham reports x = 0.65, [BT] and Bhandage reports x 
= 0.70 for Xn, 2 Mn, 2 \i ; , 2 (),. 62 The wide 
spread in critical concentrations of Mn 3+ can be ex- 
plained by the differences in JT splitting energies found 
by X-ray absorption spectroscopy on AMn 2 04 spinels by 
Noh, et al. [22] In essence, the energy drop from JT dis- 
tortion around Mn 3+ on the spinel B site is very senstive 
to changes in chemical pressure. 

We obtain the cooperative behavior of the averaged lat- 
tice using Rietveld refinement. We do not necessarily re- 
solve the distinct cation coordinations of Mg 2+ and Cu 2+ 
if they are different on the local, atomic length scale. Two 
views of the JT transition can be proposed: in the case of 
a sharp crossover, as would be implied by how Rietveld 
analysis is performed, all AO4 tetrahedra are equivalent 
whether they contain Mg 2+ or Cu 2+ , and upon increas- 
ing x they abruptly transform from ideal tetrahedra in 
the cubic spinel to flattened tetrahedra in the coopera- 
tively JT distorted spinel. In the second case, the CUO4 
tetrahedra are always locally JT distorted (even for val- 
ues of x where the spinel is cubic) but the crossover at x 
= 0.43 at room temperature represents the point where 
they cooperatively order and the JT distortions percolate 
through the long-range structure. 

Local structure via real-space analysis 

The average structure model of Mgi-zCuzC^O^ from 
Rietveld refinement indicates that the compounds exist 
as single phases, either cubic or tetragonal, apart from 
the two-phase coexistence region around Tjt- When 
modeled using a single unit cell, Mg 2+ and Cu 2+ are re- 
quired to share the same crystallographic site and their 
surroundings are necessarily identical. This model often 
inadequately describes the true structure of compounds 
where energy-lowering changes in cation coordination are 
known to persist above the average structural transition 
temperature as in perovskite manganites and cobaltites. 
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FIG. 6: (Color online) Least-squares fits to the PDF for x = 
0.20 at 300 K for the (a) cubic spinel phase with fractional- 
occupancy obtained from Rietveld refinement, and (b) an 80- 
20 fit to the end members MgC^CU and CuC^CU. 



[51 IH1 [Ml 1MJ The PDF has emerged as a key tool for mea- 
suring these local distortions that do not possess long- 
range order. Because the PDF is a weighted histogram 
of all atom-atom distances in the sample, it is sensitive 
to the distinct bond distances that are produced by dis- 
similar coordination of multiple chemical species on the 
same site. [53] We investigate whether the PDF shows 
any signature of distinct Cu 2+ and Mg 2+ coordination. 

Least-squares PDF refinements can be performed us- 
ing the average structure unit cells from Rietveld refine- 
ment as a starting point. Fits to the x = 0.20 data at 
300 K are shown in Fig.[6j Panel (a) shows the fit to a 
Rietveld-refined cubic unit cell with split 0.20/0.80 occu- 
pancy of Cu 2+ and Mg 2+ on the same crystallographic 
site. In (b), we fit using a 0.20/0.80 linear combination 
of the MgCr2 04 and CuCr204 end members with lat- 
tice parameters allowed to refine. The fit is good despite 
the use of a tetragonal unit cell to model a structure far 
above Tjt, but it does not improve on the fit using a 
single cubic unit cell. We use the value 
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to compare fits, where N is the number of points in the 
PDF from the nearest-neighbor cutoff to 20 A. The abil- 
ity to distinguish Mg and Cu is hindered by their simi- 
lar neutron scattering lengths: 5.38 and 7.71 fm, respec- 
tively. [65] For x = 0.20, least-squares PDF fits do not 
definitively prove that there are distinct (or identical) 
coordination environments for Mg 2+ and Cu 2+ . 

Resolution of distinct cation environments is aided 
when x = 0.43 due to approximately even concentra- 
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FIG. 7: (Color online) Fits to the PDF for x = 0.43 at 
300 K for (a) least-squares refinement of the two fractional- 
occupancy phases found from Rietveld refinement, (b) least- 
squares refinement using the end members MgCr2 04 and 
CuCr2 04, and (c) the fit after RMC simulation. 



tions of Cu 2+ and Mg 2+ . At 300 K, Rietveld refinement 
found coexistence of the cubic and tetragonal phases. 
Neither phase alone can be used to produce a satisfac- 
tory fit to the PDF. A two-phase fit using the Rietveld 
refined cells is shown in Fig. [7^ a), and agrees quite well 
with the data. As with the x — 0.20 sample, we also fit 
the data to a combination of the end members MgCr204 
and CuCr 2 4 in Fig.^b). A gain, the fractional occu- 
pancy Rietveld result produces a better fit than the end 
members. The least-squares fits indicate that the average 
structures produce excellent representations of the local 
structures, but they do not definitively show whether the 
Mg 2+ and Cu 2+ coordination environments are distinct 
or similar. 

Least-squares fits to the PDF are required to spec- 
ify how many distinct AO4 environments to allow, much 
like in a Rietveld refinement. The Fdim and IAi/amd 
unit cells provide only one AO4 environment per phase. 
There is no way to define cation-dependent coordination 
without manually building a lower-symmetry unit cell. 
In order to investigate the AO4 environment directly, we 
remove the symmetry constraints of least-squares PDF 
analysis and utilize large-box modeling via reverse Monte 
Carlo (RMC) simulations. This method has proved to 
be useful for investigating atomic structure on the local 
level, especially in cases where long-range periodicity is 
not present, such as SrSn0 3 ,[2S] Bi 2 Ti 2 07,[2I] and /3- 
cristobalite. [BB] 

The RMC supercell for an x = 0.20 sample at 15 K 
is inspected by folding each of the unit cells back into a 
single box, shown in Fig.JHJ which reveals how each crys- 
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FIG. 8: (Color online) The final 10 x 10 x 7 RMC supercell 
for an x = 0.20 simulation is folded into a single unit cell 
to obtain point clouds at each crystallographic site. Atoms 
shown are Cu (red), Mg (dark blue), Cr (light blue), and O 
(orange). 
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FIG. 9: (Color online) Bond valence sums for Cu and Mg 
extracted from Mgi_ a; Cu :l ;Cr204 RMC supercells x = 0.20 
after fits to data at T = 300 K (a) and 15 K (b). Aside from a 
slight broadening at 300 K, the distributions are comparable. 



tallographic site is decorated with atoms. The supercell 
contains 560 Cu and 2240 Mg atoms that are randomly 
arranged. Because there are a large number of distinct 
Cu and Mg atoms, statistical analysis can be used to in- 
vestigate whether there is any evidence for the local AO4 
distortion to depend on the central cation. Bond valence 
sum histograms shown in Fig. [9] show that both the A 
cations have valences peaked around the expected value 
of A 2+ , and the B site shows only Cr 3+ , which implies 
that our supercells contain chemically reasonable bond 
lengths. 

No cation dependence of AO4 bond distances is ob- 
vious from partial radial distributions gcu-o( r ) an d 
9M g -o{ r )- This is to be expected because Mg 2+ and 
Cu 2+ have the same ionic radius when tetrahedrally co- 
ordinated. Neither is there any apparent distinction be- 
tween Cu0 4 and Mg04 tetrahedra based upon O—A— 
O bond angles, in contrast to our previous study on 
CuMn 2 4 .[23] 

We use the continuous symmetry measure (CSM) tech- 
nique to gauge the tendency for JT distortion of AO4 
tetrahedra. The CSM technique provides a distance mea- 
sure (DM) of a given tetrahedron that indicates its de- 
viation from ideality. [301 131] A perfect tetrahedron has 
DM = 0, and any distortion increases the value of DM. In 
the end member compounds, Mg0 4 has a DM = while 
Cu0 4 has a DM = 0.0076. Having thousands of distinct 
tetrahedra in the RMC supercell affords the opportunity 
to produce a histogram of DM (Fig.fToj) for all tetrahe- 
dra depending on the central cation. Each panel contains 
eight lines: four for each cation, resulting from four in- 
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FIG. 10: (Color online) CSM histograms from RMC simu- 
lations for three Cu contents and three temperatures. Each 
panel contains 8 distributions: one Cu and one Mg for four 
independent simulations. Broadening across the T = 300 K 
series indicates increasing tetrahedral distortion with x. Ther- 
mal effects lead to sharpening and shifting toward smaller DM 
when the x = 0.20 sample is cooled to 150 K. This tendency 
toward ideal tetrahedra is disrupted by long-range tetragonal 
distortion upon further cooling, so the T — 15 K histogram 
is broad and shifted to higher DM. 
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FIG. 11: (Color online) Cumulative CSM distributions are 
similar enough from run to run that most lines overlay each 
other. Only the x = 0.43 run has a clear distinction between 
CUO4 and Mg04 tetrahedron shape, with Mg04 distinctly 
more ideal than CuCU. 
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FIG. 12: (Color online) CSM distributions of AA4. tetrahedra 
(connected by long red bonds in inset) as a function of number 
of neighbors that are JT-active Cu 2+ . Distinct distributions 
are shown for CUA4 (circles) and Mg J 44 (diamonds). Linear 
guides to the eye for each distribution show increasing tetra- 
hedral distortion with number of Cu neighbors, and overall 
more distorted tetrahedra when tetrahedra are Cu-centered. 



dependent RMC simulations. The overall shape of each 
histogram describes the average distortion (peak center) 
and the tightness of the DM distribution (peak width) at 
each value of x and T. No tetrahedra are present with ex- 
actly DM = because stochastic RMC simulations leave 
no atomic positions untouched-even a compound with 
ideal tetrahedra would have the shapes subtly distorted. 

At 300 K, the spread of distortions increases for x = 
0.43 because the average structure becomes a mixture of 
cubic and tetragonal phases. In the series where x = 
0.20, the peak sharpens upon cooling to 150 K, which 
we attribute to a reduction of thermal vibrations. It also 
moves to smaller DM values, indicating a tendency to- 
ward a more ideal MO4 environment. This trend would 
continue to low temperature in the absence of long-range 
JT distortion. Instead, the peak broadens and shifts to 
higher DM at 15 K. We attribute both of these effects 
to the tetragaonal phase transition. These histograms 
provide a view of the average tetrahedral shape, but no 
distinction between Mg and Cu is apparent. We find that 
cumulative distributions offer a clearer picture of this de- 
pendence. 

The cumulative distributions in Fig. [IT] do not show 
separation between the Cu and Mg curves in the x = 
0.10 or 0.20 samples. However, the x — 0.43 sample 
shows a clear distinction, with Mg tetrahedra possessing 
DM that are closer to zero (more ideal) than Cu. This 
is clear evidence for the tendency of Cu to undergo JT 
distortion while Mg remains more symmetric. 



Overlapping DM curves for x — 0.10 and 0.20 do not 
preclude the possibility of distinct CUO4 and Mg04 tetra- 
hedra in those samples. There may be insufficient res- 
olution in the PDF to distinguish the cations due to 
their similar neutron scattering cross sections and the 
relatively low concentration of Cu 2+ . It remains un- 
clear whether distinct cation environments are only seen 
around x = 0.5 (due to approximately even cation con- 
centrations) and if the distinction would disappear as 
x — >• 1. Temperature dependence of the CSM when x 
= 0.43 may also provide some insight into the dynamics 
of these distortions. Substituting a JT-inactive A site 
cation with a different neutron scattering cross length, 
such as Mn 2+ or Co 2+ , may aid contrast with Cu 2+ . 

Short-range interactions between JT-active cations are 
evident in the A-A correlations in Fig. [12] Each A cation 
has four A nearest neighbors 3.59 A away, arranged in a 



diamond lattice (inset in Fig. 12 1. The four nearest neigh- 



bors create a large tetrahedron around the central cation. 
The tetrahedral DM for these AA4 are plotted for each 
central cation (Mg or Cu) as a function of the number 
of Cu nearest neighbors, from to 4. Multiple RMC 
simulations show an upward trend indicating more dis- 
tortion as the number of nearby Cu increases. The sep- 
aration between the two lines implies that Cu-centered 
Cu^4 have higher tetrahedral DM (more distortion) than 
Mg-centered MgAt. These short-range correlations are 
indicative of strain coupling between adjacent A-site JT 
distortions. 
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FIG. 13: (Color online) Magnetic hysteresis of 
Mgi-^CuzC^CU with x = 0.43 at T = 5 K after zero- 
field cooling (ZFC) and field cooling (FC) with H FC = 5 T. 
The shift in the —H direction when M = is the exchange 
bias field He- 



Magnetic properties 

MgCr204 and CuC^C^ have markedly different mag- 
netic behavior due to the addition of unpaired spins in 
Cu 2+ and the accompanying JT distortion. CuCr204 is 
a hard ferrimagnet with Tc = 135 K,[67 while MgC^C^ 
undergoes complex antiferromagnetic ordering below Tjy 
= 16 K. 53, 68J Macroscopic composites of a ferromag- 
net and antiferromagnet would result in a traditional ex- 
change biased material, with an enhanced coercive field 
He and an exchange bias field He, manifested as a shift 
of the hysteresis loop in the —H direction. [55H7T] 

The Mgi^zCuzC^Cu solid solution is a mixture on 
the atomic level, but nevertheless exhibits the magnetic 
hallmarks of an exchange biased system. The hystere- 
sis loop of the x — 0.43 sample in Fig. [13] at T = 5 K 
has H c = 8.4 kOe. Cooling with a field H FC = 50 kOe 
broadens the hysteresis loop and shifts it in the —H direc- 
tion by He = 0.44 kOe. This shift signifies the preference 
for the ferrimagnet to align along the field-cooling direc- 
tion. He and He decrease with the Cu 2+ concentration 
x , as does the onset of magnetic ordering. These trends 
In all cases, field-cooling increases 



are shown in Fig. 14 



He and results in the appearance of a significant He- 
We find He = in all samples after zero-field cooling. 

Exchange bias is traditionally manifested by pinning at 
the interface between a ferro- or ferrimagnet and antifer- 
romagnet after field cooling when Tc > T^. It can also 
arise from intrinsic disorder in a single-phase system, or 
the presence of disordered spins from either a spin glass 
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FIG. 14: (Color online) Coercive fields He (left) and exchange 
bias fields He (right) for Mgi_ x Cu :c Cr204 samples with x = 
0.43, 0.20, and 0.10 after ZFC and FC with H FC = 5 T. In 
all cases, ZFC results in a smaller He and gives He = 0. 



or uncompensated surface spins on small particles (which 
behave in a glassy manner themselves). [73] In the 
Mgi-zCuzC^O^ system, the magnetic behavior cannot 
be fully described by traditional AFM-FM interplay be- 
tween MgCr204 and CuCr204 because the onset of He 
when x = 0.43 in Fig.[T4]occurs above T = 20 K, which is 
above the Neel temperature of MgCr 2 04. Therefore the 
disordered solid solution contains some regions of het- 
erogeneity which may behave as glassy moments or as 
some intermediate AFM phase. Either of these cases can 
produce exchange bias. [HI [74] 

The correlations between A cations seen in Fig. [^indi- 
cate that the local structure varies with the composition 
of nearby atoms. Magnetic interactions, in turn, will 
be affected by these local distortions. A more detailed 
investigation of magnetic behavior, as was performed on 
Zn :E Mn3_ :E 04, j55j may help elucidate how the competing 
structural and magnetic interactions produce exchange 
bias. The phenomenon of JT-active cation clustering has 
been investigated by dilatometry, [75] but its effects on 
magnetism have not been explored in detail. 



CONCLUSIONS 

We find that the solid solution Mg^^-Cu^C^C^ has a 
two-phase coexistence of cubic and tetragonal phases at 
room temperature for 0.43 < x < 0.47. Tetragonality is 
induced by increasing JT activity in tetrahcdrally coordi- 
nated Cu 2+ . The average structure descriptions from Ri- 
etveld refinement provide an adequate description of the 
structures across the range of x and temperature. This is 
corroborated by the magnetic behavior, which indicates 
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a disordered atomic mixture. The x — 0.20 sample is 
cooperatively JT distorted with orbital ordering at 15 K 
despite 82% of all Cu 2+ having zero or only one Cu 2+ 
neighbors. 

Least-squares PDF refinements achieve good fits us- 
ing the models from Rietveld refinement, implying that 
it might be difficult to improve on models where CuC>4 
and MgC>4 are equivalent. Still, bond valence calcula- 
tions show that RMC simulations produce better fits to 
the data while retaining chemically reasonable bond dis- 
tances. The AO4 tetrahedral distortion increases with x 
as judged by CSM. The absence of distinct CUO4 local 
distortions in the low- a; CSM histograms does not pro- 
hibit their existence. When x — 0.43, tetrahedral DM 
indicate distinct coordination of Cu 2+ versus Mg 2+ . Fur- 
thermore, cation-cation interactions probed by AA4 CSM 
indicate that local clustering of Cu leads to increased JT 
distortion. This technique resolves cation-dependent JT 
distortions (even when they are incoherent) in materials 
where contrast between cations exists. The presence of 
magnetic exchange bias implies that short-range struc- 
tural details are influencing the magnetic interactions, 
and more complex magnetic characterization may help 
describe these interactions. 
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